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I.  INTRODUCTION 


In  Reference  1  we  presented  a  derivation  of  both  the 
electric  field  integro-dif ferential  equation  (EFIDE)  and  the 
magnetic  field  integral  equation  (MFIE)  for  the  determination 
of  the  current  density  induced  on  the  surface  of  metallic 
(perfectly  conducting)  structures  situated  in  a  vacuum  half¬ 
space  above  a  lossy  half-space.  We  showed,  by  comparison  to 
experiments  performed  at  the  University  of  Michigan,  the  quality 
of  the  predictions  that  could  be  obtained  by  using  the  patch 
zoning  method  to  solve  the  MFIE  for  two  limiting  cases  of  the 
lossy  half-space:  free-space  and  perfectly  conducting.  The 
code  developed  for  this  purpose  is  capable  of  obtaining  surface 
current  predictions  for  an  aircraft  model  consisting  of  inter¬ 
secting  elliptical  cylinders  (Figure  1) ,  optionally  in  the 
presence  of  a  perfectly  conducting  ground  plane.  Reference  1 
presents  comparisons  of  predictions  versus  experiments  for  the 
full  aircraft  model  situated  in  free-space  as  well  as  a  cylinder 
above  a  perfectly  conducting  half-space. 

In  this  report  we  present  the  result  of  our  efforts  to  enhance 
our  computer  code  (NEC-2A)  by  including  the  capability  to  treat 
the  same  aircraft  model  in  the  presence  of  a  lossy  half-space. 
Before  embarking  on  the  extensive  programming  effort  required 
to  solve  the  MFIE  for  the  two-layered  problem,  we  checked  the 
validity  of  our  equations  by  rederiving  them  in  an  alternate 
manner.  We  present  this  derivation  here,  as  well  as  its  general¬ 
ization  to  multilayered  horizontally  or  radially  stratified 
media.  We  show  that  the  electric  dipole-excited  electric  or 
magnetic  fields  are  the  fundamental  quantities  required  to 
determine  the  EFIDE  and  the  MFIE  respectively.  We  utilize 
these  dipole-excited  quantities  to  construct  Green's  dyadics, 
and  the  procedure  for  doing  this  is  presented  in  detail.  The 
determination  of  these  dipole-excited  fields  is  basically  as 


GROUND  PLANE 


Figure  l.  Aircraft  Model 


tractable  for  the  multilayered  problem  as  it  is  for  the  two¬ 
layered  problem.  Following  a  procedure  presented  in  the  book 
by  Felsen  and  Marcuvitz  (Reference  2) ,  one  can  determine  all 
the  necessary  dipole-excited  fields  from  scalar  potentials. 
Utilizing  these  potentials  and  satisfying  boundary  conditions 
leads  to  a  system  of  linear  algebraic  equations  for  a  finite 
number  of  unknown  constants  and  this  number  depends  linearly 
on  the  number  of  layers  used  to  describe  the  medium.  If  we 
intend  to  develop  a  computer  code  for  the  use  of  the  dipole 
solutions,  e.g.,  the  MFIE,  then  the  code  can  be  written  to 
solve  the  finite  system  of  equations  as  part  of  its  execution. 

To  make  the  existence  of  the  dipole-excited  solutions  have 
more  than  just  a  potential  role  for  future  work,  we  present 
the  explicit  use  of  these  dipole  solutions  (Reference  3)  for 
the  horizontally  stratified  two-layered  medium  to  obtain  exactly 
the  same  EFIDE  and  MFIE  derived  in  Reference  1  and  thus  complete 
the  desired  alternate  derivation.  Further  use  of  our  general  . 
formulation  can  be  made  by  using  both  the  analysis  and  computer 
code  described  in  Reference  4  for  the  horizontally  stratified 
three-layered  medium.  Additional  material  on  dipole-excited 
fields  in  a  multilayered  medium  can  be  found  in  Reference  5. 

We  exercised  certain  switches  in  the  computer  code  developed 
for  the  two-layered  problem  to  compare  our  code's  capabilities 
to  two  other  sets  of  data  for  the  response  of  a  perfectly  con¬ 
ducting,  finite  length  cylinder  in  the  proximity  of  planar 
lossy  material.  One  set  of  data  is  the  current  density  measured 
on  the  cylinder  when  it  was  placed  in  front  of  lossy  material 
within  the  anechoic  chamber  at  the  University  of  Michigan 
(Reference  6) .  The  other  set  of  data  furnished  to  us  by  G.  J. 
Burke  was  calculated  as  a  special  case  of  a  Lawrence  Livermore 
Laboratory  code's  capability  which  corresponds  to  a  special 
case  of  our  code's  capability.  A  more  detailed  description  of 
that  code  can  be  found  in  Reference  7.  This  data  was  the  plane 
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wave  induced  bulk  current  on  a  finite  cylinder  above  and 
parallel  to  a  lossy  half-space.  The  results  of  both  the 
Michigan  and  Livermore  comparisons  as  well  as  other  inter 
mediate  tests  are  presented  in  this  report. 
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II.  SURFACE  INTEGRAL  REPRESENTATIONS  FOR  SCATTERING 
BY  AN  OBJECT  EMBEDDED  IN  A  STRATIFIED  MEDIUM 

The  presentation  of  this  analysis  is  facilitated  by  refer¬ 
ring  to  Figure  2.  In  that  figure  we  divide  all  of  space  into 
N  distinct  volumes  and  label  the  volume  containing  the  object 

with  the  index  N  .  Each  volume  is  characterized  by  different 
s 

electrical  parameters  so  that  Maxwell's  equations  appropriate 
for  a  particular  volume  V£  are  written  as 

V  x  E£a(r)  *  iu)U^H^a  (r)  (1) 

V  x  H£a(r)  =  -ia)e£E£o(r)  +  J£a(r)  (2) 

where  Z-1,2 , . . . ,N  and  for  the  case  that  will  be  developed  in 
most  detail 


eZ  ~  eZR  +  ieZ I  (3) 

The  rigid  source  terms  satisfy  the  condition 

J£ot(r)  ji  0  only  if  reV£ocv£  (4) 

and  in  this  section  a=a  will  correspond  to  mathematically 
introduced  delta  function  sources  and  a»b  will  correspond  to 
the  physical  distribution  of  sources.  The  analysis  that  will 
be  presented  has  applicability  for  more  general  situations 
than  the  stratified  model  of  all  space  depicted  in  Figure  2; 
however,  that  situation  includes  all  cases  of  current  interest. 

When  the  sources  are  delta  function  sources  located  at  a 
point  r a,  then  we  introduce  the  common  notation 
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1 


=  St<£.£ a) 

=  £a{£»£o) 


(5a) 

(5b) 


Whenever  the  sources  are  delta  functions,  we  will  consider  the 
scattering  body  to  be  absent.  We  now  use  the  vector  identity 

7 * ( AxB )  =  B* (VxA)  -  A* (VxR) 


to  obtain 


7  *  (— f.bx-£a^  ”  a7*-^b 

-  (  *ti**h&ti**Bu)  <6) 

The  quantities  E^a  and  H^a  satisfy  Equations  1  and  2  for  a=a, 
and  similarly  E£b  and  satisfy  Equations  1  and  2  for  a*b. 
Substituting  these  versions  of  Equations  1  and  2  into  the 
right-hand  side  of  Equation  6  and  then  cancelling  terms,  we 
obtain 


-lb*  -£a"-£a*-£b 


(7) 


D„  =  E.  axH_. -E..  xH. 
-*  -£a  -£b  -£b  -£a 


(8) 


Integrating  both  sides  of  Equation  7  over  and  employing 
the  divergence  theorem,  we  obtain 


Is£,£-1  +  Zs£,£+1  +-  5£,N  Zs  ZV£ba  "  ZV£ab 

s 


(9) 


where 


Zs£,k 


(10) 


13 


(11) 


i.  =  f 


/• 

0  ~  D  =  J 


V£a6 


P  -N. 


E.  .J  dV 

-&a  —IQ 


where  is  the  surface  interface  between  volumes  V£  and  V^; 

n^^  is  the  normal  to  ^  directed  outward  from  into  V^; 

S  is  the  surface  of  the  perfectly  conducting  object  which  is 

embedded  in  V„  ;  ft  is  the  inward  normal  to  S.  and  6„  „  is 
Ng  P  ^,NS 

the  Kronecker  delta  function.  In  Appendix  A  it  is  argued 

that  despite  the  absence  of  physical  surfaces  for  I  n  and 

IgN  N+1/  both  integrals  have  zero  value. 

We  now  focus  attention  on  Equation  9  for  the  situation 

where  £  =  Ng;  the  only  J^a  source  is  in  VN  and  has  the  form 

s 


ft  6  (r-r^ 


For  this  case.  Equation  9  becomes 


where 


FN  ab  +  Xs  "  a  '  -N  b  *£a* 
s  s 


FNgab  =  IsNg,Ns-l  +  IsNg,Ns+l  +  IVNgab 


In  Appendix  A,  it  will  be  shown  that  FN  afc  is  independent  of 

s 

the  presence  (or  absence)  of  the  embedded  perfect  conductor 

for  a  general  distribution  of  J^'  s,  and  this  proof  will 

utilize  Equation  9  for  i  not  being  restricted  to  N  .  Using 

the  fact  that  FN  ab  does  not  depend  on  the  scattering  object 

s 

we  will  give  this  quantity  explicit  meaning.  We  return  to 
Equations  13  and  14  and  consider  that  we  sequentially  intro¬ 
duce  the  following  three  sources  all  located  at  r_ 


£  a  -  Mr-^) 


i=l #2,3 


(16) 


with  ( a^,a2,a forming  an  orthonormal  set,  and  for  each  i 
we  obtain,  in  the  same  manner  that  Equation  14  followed  from 
Equation  13,  the  relation 


FNgab  +  Is  ~  ai *  ^Ngab 


(17a) 


where 


ab  =  b(^aJ 
s  s 


(17b) 


Utilizing  the  independence  of  the  unit  vectors,  we  rewrite 
this  set  of  equations  as  a  single  vector  equation. 


where 


^Ngab  +  ^s  *  al(al*^sab) 

+  a2^2'^Nsab| 

+  a3  (a3*  ^Ngabj  (18) 


F  =  £ 

i-Ngab  "  l  Ngab 


a  A  aa  2  A  aa  3 

2FN _ab+  a3FN  ab 
s  s 


(19) 


and 


Is 


-  Mi  +  M ?  +  a,i3 

is  2  s  3  s 


(20) 
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The  right-hand  side  of  Equation  18  can  be  interpreted  by 
noting  that  the  identity  dyadic  I  is  expressible  in  terms 
of  the  set  (a ^,a2>a3)  as 


I  *  aia1  +  a2a2  +  a3a3 


so  that 


ll(ai*?Nsab)  +  a2(a2*^Nsabj  +  S^'^abJ 


i*^N  ab 


‘=N  ab  -N  ab 
s  s 


and  Equation  18  becomes 


^N_ab  +  ~  ab 

s  s 


Using  the  fact  that  FN  ab  is  independent  of  S,  we  consider 

s 

the  special  case  of  Equation  23  where  S  vanishes  and  in  turn 
I  vanishes,  so  that 


ab  “  ^TN  b  -a*  ^TN  ab 
s  s  s 


with  .  being  the  total  electric  field  at  r  excited  by 

* tn_ ao  -  a 

s 

a  general  distribution  of  sources  ETN  includes  all 

s 

multiple  scattering  effects  due  to  the  stratified  medium  but 
does  not  include  any  scatterng  due  to  the  embedded  perfect 
conductor.  We  now  write  Equation  23  as 


^Nab<-a)  "  5TNgb(-a)  +  -s 
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We  now  focus  our  attention  on  1^  and  express  Equation  20  in 
more  detail 


I 

— s 


'•^hds 


(26) 


with 


dJ  =  Ei  a  x  H..  . 
— N  — N  a  — N  b 

s  s  s 


-  'Sh  b  x  £  a 
s  s 


(27) 


Forming  n  *D^  and  using  standard  triple  scalar  product  mani- 
p  — «s 

pulations,  we  obtain 


n  ‘D1 
p  -N 


s 


•  ( n  xH  ‘E^l  -/n  xE„  .  \ 

\  P  “Nsb/  “V  \  p  ~®sh) 


(28) 


The  boundary  conditions  require  that  the  tangential  electric 
field  excited  by  the  J^'s  on  the  embedded  conductor  vanishes 
so  the  second  term  in  Equation  28  vanishes.  We  note  in  passing, 
that  we  cannot  apply  the  same  argument  to  eliminate  the  first 
term  since  the  scattering  body  is  absent  in  the  "a"  system. 
Equation  28  is  now  written  as 


n 


•D1 
P  -N 


s 


-5Ssa 


(29) 


In  order  to  obtain  this  equation,  we  used  the  fact  that  the 
surface  current  on  S  excited  by  the  s  is  given  by 

J(r)  -  *„  x  b  (30) 

s 

with  ftp  *  -fiQ  and  nQ  is  the  outward  normal  to  S .  Combining 
Equations  25,  26,  and  29,  we  obtain 

-&mb(£a>  SN  <£•£»><“  ‘W 

s  s  s  s 
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where 


(— '-a 


(32) 


Equation  31  has  many  of  the  properties  we  were  seeking  in  that 
it  is  a  representation  of  the  fields  excited  by  general 
sources  in  a  stratified  medium  having  an  embedded  conductor; 
the  presence  of  the  conductor  is  accounted  for  by  the  sur¬ 
face  integral.  Furthermore,  this  surface  integral  contains  a 
Green's  dyadic  that  can  be  explicitly  formed  in  the  described 
manner,  i.e.  using  the  three  electric  field  representations 
generated  by  three  orthonormal  dipole  sources.  Despite  these 
features  of  Equation  31,  it  is  not  the  final  form  desired 
because  the  Green's  dyadic  appears  to  be  the  right  of  J(r) , 
and  this  form  does  not  readily  allow  us  to  take  the  curl  of 
both  sides  of  the  equation.  The  significance  is  that  this  is 
the  procedure  used  to  obtain  a  surface  integral  representa¬ 
tion  for  H£ (r) . 

Returning  to  Equation  31,  we  use  the  following  property 
valid  for  any  dyadic  D  and  vector  V 

V-D  *  D-V  (33) 


where  D  is  the  transpose  of  D.  The  surface  integral  in 
Equation  31  is  now  written  as 


J  (£'£a)ds  *  J  iN  (£*£3) •  J(r)dS  (34) 


s  s 


In  Appendix  B  we  will  prove  that 
gj (E'Sa>  =  §k<ra,r) 


^  eVj  '  -aeVk  (35) 


Combining  Equations  31,  34,  and  35,  we  obtain 


18 


(36) 


— N  b (— a>  =  b^a*  +J  &  <£a'£)'2(£,ds 
s  s  s  s 

which  is  our  desired  representation  for  the  electric  field. 
Next  we  introduce  a  notation  change  in  order  to  make  this 
equation  conform  to  standard  notation  for  integration  and 
observation  variables;  i.e.,  r  is  typically  used  for  the 
observation  variable  and  r'  for  the  integration.  The  distinc- 
tion  between  r  and  r*  is  more  than  just  notation  in  that  r 
must  be  interpretable  as  ranging  over  a  three-dimensional 
volume  while  r’  is  restricted  to  range  over  the  two- 
dimensional  surface  of  the  scattering  object.  Denoting  the 
volume  variable  in  Equation  36  as  r  and  the  surface  variable 
as  r'  and  dropping  the  explicit  "b"  notational  dependence, 
we  rewrite  that  equation  as 

K  (r)  -  E  (r)  +/  Gn  (r  ,r ' )  *  J(£' )  dS '  (37) 

-“s  s  s  s 

Before  we  elaborate  on  the  meaning  of  this  representation,  we 
claim  that  the  representation  is  valid  for  general  values  of 
1;  and  for  general  i,  we  have 

Efc(r)  -  5r*<£)  +  f  (£»£')  *J(r)dS'  (38) 

s 

The  essential  details  of  how  Equation  38  is  derived  for  £/Ng 

is  presented  in  Appendix  C.  The  reason  that  we  focus  on  the 

case  £=N  is  that  it  is  this  case  that  leads  to  the  explicit 
s 

equation  that  we  numerically  solve  to  determine  J(r) . 

A  discussion  of  the  meaning  of  Equation  38  requires  that 
we  return  to  the  definition  of  given  by  Equation  32.  With¬ 
out  restricting  1,  we  use  the  arguments  leading  to  that 
equation  to  express  (r,r')  as 


L9 


3 

20<E'E,)=£  Ej(r,r')a.  (39) 

l  i«l  *  -  * 

with  E^(r,r’)  satisfying  the  boundary  conditions  described 
in  Appendix  A  as  well  as  the  equations 

V  x  E*(r,r')  =  iwu£Hj(£,r’) 

V  x  H^(r,r’)  =  -iaje^E^  (r,r 1 )  +  a±  6(r-r') 

In  these  equations  r  is  the  variable  ranging  over  a  three- 
dimensional  space,  V  ,  and  the  curl  appearing  in  these  equa¬ 
tions  operates  with  respect  to  this  variable.  The  r'  pertains 
to  the  location  of  the  dipole  source,  and  it  is  chosen  to  lie 
in  the  same  volume  as  the  embedded  scatterer.  It  is  oriented 
in  the  sL  direction  with  the  choice  of  these  orientations 
being  such  that  (a1,a2»a3)  form  an  orthonormal  set.  To  com¬ 
plete  the  description  of  Equation  38,  we  simply  state  that 
ET^(r)  is  the  total  electric  field  at  r,  which  includes  all 
multiple  scattering  due  to  the  stratification  of  the  medium, 
but  considers  the  embedded  conductor  to  be  removed. 

Next  we  obtain  a  representation  for  the  magnetic  field  in 

taking  the  curl  with  respect 
of  Equation  38  to  obtain 

V  x  G^(r,r') •J(r')dS  (42) 

case  o  *  b  and  considering  the 
we  would  be  able  to  conclude  that 

iuiM^H^  (r)  (43) 


terms  of  the  surface  integral  by 
to  the  r  variables  of  both  sides 


V  x  EA(r)  =  V  x  E^U)  +  / 

s 


Using  Equation  1  for  the  special 
scattering  object  to  be  present, 


V  x  E£(r) 


(40) 

(41) 
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or  considering  the  object  to  be  absent,  we  would  be  able  to 
conclude  that 


7  x  58  ia)U«,5TJi<£) 


(44) 


Now  defining 


' )  s  7  x  £*,<£'£')  (45) 

we  can  rewrite  Equation  42  as 

(r)  =  (£)  +  y*KA(£,£^ ')  *J(r ')dS'  (46) 

s 

We  obtain  our  final  expression  for  by  combining  Equations 
39,  40,  and  44  which  results  in 

3 

K^r’)  =  £  Hj(r,r’)  a±  (47) 

i-1 

Equations  38  and  46  are  the  desired  surface  integral  repre¬ 
sentations  and  Equations  39,  40,  41,  and  47  are  the  equations 
that  give  these  representations  utility.  In  the  next  section 
we  will  given  Equations  39  and  40  explicit  meaning  for  a  two¬ 
layered  medium,  and  in  the  following  section  we  will  obtain 
the  EFIDE  and  the  MFIE  by  discussing  the  limit  of  Equations 
38  and  46  as  r  approaches  the  surface. 
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III.  EXPLICIT  GREEN'S  DYADIC  REPRESENTATIONS  FOR 
THE  TWO-LAYERED  CONFIGURATION  OF  INTEREST 


In  this  section  we  will  obtain  explicit  representations 
for  (£,£*)  and  (£»£')  using  the  prescription  summarized 
by  Equations  39,  40,  41,  and  47  with  the  boundary  conditions 
described  in  Appendix  A  also  being  satisfied.  These  repre¬ 
sentations  will  be  found  for  the  case  where  we  have  a  two¬ 
layered  horizontally  stratified  medium  with  one  layer  being 
a  lossy  half-space,  £=1,  and  the  other  layer  being  a  vacuum 
half-space,  1-2.  Our  scattering  object  is  in  V^»  and  we  will 
focus  our  attention  on  obtaining  the  explicit  representations 
for  S2  (r ,£*)  and  g2  (£»£')•  We  present  these  represen¬ 

tations  in  a  form  where  it  is  clear  that  they  are  identical 
to  representations  of  the  same  dyadic s  that  were  obtained  by 
a  different  formalism  {Reference  1) .  We  obtain  our  explicit 
representations  by  utilizing  representations  contained  in  the 
book  by  Banos  (Reference  3}  for  the  electric  fields  and  the 
magnetic  fields  excited  in  the  vacuum  half-space  above  a 
lossy  half-space  by  horizontal  and  vertical  electric  dipoles 
located  in  the  vacuum  region. 

% 

In  order  to  use  Banos's  explicit  results,  it  is  necessary 
to  employ  a  systematic  bookkeeping  procedure  to  compensate 

f\j 

for  the  choice  of  coordinate  system  chosen  by  Banos,  which  was 
convenient  for  his  purposes  but  not  ours.  Specifically,  Banos 
gives  explicit  solutions,  that  we  will  denote  E02(rD,rp  and 
— B2  (— d'£'  D^  wh*-ch  satisfy  the  equations 


7D  X  5b2^^D'— *  *^2— B2 


(48) 


VD  *  ^B2 (£d'£dJ 


"iwe2^B2  +  P*i*(£D“S^>  (49) 
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and  it  is  necessary  to  augment  the  solutions  to  these  equations 

for  two  reasons  in  order  to  obtain  the  £2 (£»£')  and  H^^r') 

required  by  Equations  39  and  47.  One  reason  is  that  Banos  chooses 

a  coordinate  system  in  which  r^  is  allowed  to  vary  only  over 

one  linear  dimension  rather  than  the  necessary  three,  and  the 

other  is  that  he  presents  explicit  solutions  for  only  two 

orthonormal  vectors  =  a  ,  and  a,  =  a  .  Despite  these 

i  x  J  z 

limitations,  the  necessary  solutions  can  readily  be  obtained 
from  E^'3  and 

We  now  elaborate  on  the  coordinate  system  used  by  Banos. 

The  origin  of  that  coordinate  system  lies  on  the  interface, 
and  the  positive  z-axis  points  into  V_.  The  subscript  D  is 

'Vi  ^ 

not  used  by  Banos  but  is  introduced  by  us  in  order  to  dis¬ 
tinguish  between  the  variables  used  in  the  previous  section 
from  his  variables.  We  use  his  result  by  obtaining 

— B2  D'~ and  — B2^^D'— D*  an<*  then  utilizing  the  following 

relationship 

4(£'E'>  -  |y 

S2<£.£->  -  [rD(r,r'l  ,  r'(r,r')] 


(50) 

(51) 


In  order  to  use  Equations  50  and  51,  it  is  necessary  to 
determine  the  functional  relationships 


(52) 


and 


.1 

=-D 


r-  =  £3  (r,r*  ) 


(53) 
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The  description  of  how  the  coordinate  functional  relation¬ 
ships  are  obtained  is  facilitated  by  referring  to  Figure  3. 


This  figure  is  the  view  of  the  z  =  0  plane  (the  interface) 
for  all  coordinate  descriptions  of  interest.  Because  of 
this  choice  of  coordinates,  it  follows  that 


zD  »  z  (54) 

*D'  =  z‘  (55) 
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and  we  note  that  Banos  uses  the  notation  z^  =  h.  The  signifi¬ 
cance  of  choosing  the  coordinate  system  having  the  origin  "0" 

.  'V* 

in  addition  to  the  Banos  coordinate  system  having  the  origin  ”D' 
is  that  the  location  of  the  dipole  source  is  fixed  in  the  HDM 
coordinate  system  at  the  origin  but  is  free  to  vary  in  the  "0" 
coordinate  system.  Also,  its  position  is  denoted  £’  in  accord 
with  the  meaning  of  the  primed  coordinates  described  in  the 
previous  section.  Still  in  accord  with  the  meaning  ascribed 
to  the  variables  discussed  in  that  section,  £  identifies  the 
location  of  the  observation  point  in  the  "0"  coordinate  system. 
From  the  construction  of  the  two  coordinate  systems,  we  have 
the  relationship 


£o  =  £-£  (56) 

and  the  fact  that  p_  is  the  difference  between  the  two  coordi- 
nates  of  interest  is  the  reason  for  designating  Banos  coordi- 
nates  with  the  subscript  “D".  It  should  be  noted  that  Banos 
denotes  pQ  =  jp^J  as  r  in  his  book,  and  this  should  not  be 
confused  with  jrj  where  we  mean  r  =  £  +  z We  now  have 
determined  the  required  coordinate  functional  relationships 
which  can  be  summarized  as 

£d  =  £-£  +  zSz  (57) 

r^  *  z’  az  (58) 

2  2 

We  will  now  determine  E_-(r  ,rl)  and  H__(r  ,rM  in  terms 
of  the  explicit  expressions  that  Banos  gives  for  (f^r^) 
and  3b2 *  T^e  derivation  simplifies  by  noting  that  the 
spatial  dependence  of  the  latter  quantities  can  be  expressed 
in  terms  of  the  three  scalars  pD»  z+z’,  and  $D.  Let  us  con¬ 
sider  the  situation  described  by  the  argument  (pD,z+z*  ,$D+ir/2) 
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and  the  fields  excited  by  an  electric  dipole  located  at  the 
°rigin  of  the  "D"  coordinate  system  that  is  directed  along 
=  ^2*  The  fields  generated  by  the  dipole  are 

-B2(PD'Z+Z'  ^D+1t/2>=  EpD2{pD'Z+zf  '4,D+ir/2)apD(4’D+1T/2) 

+  E^B2(pD'Z+Z'  '  Vir/2)%D(V1f/2) 

+  EzB2(pD'Z+Z,'VV2)&z  (59) 

The  geometry  of  the  physical  situation  is  such  that 

EpB2(PD'Z+Z’'V’,/2)  ’  EJb2(PD'Z+2''V  ,60a) 

E»B2(PD'2+2'  'V721  *  E*B2(VZ+Z''V  <60b) 

EzB2(pD'2+2‘'V”-/2)  *  EzB2(PD'Z+2''V  (60c1 

and  it  is  worth  noting  that  despite  these  relationships 

£b2<PD'2+z''V,'/2)  *  Sb2!pD'2+z''V 

because  { d>D+-rr/2)  /  Sp(<}>D),  and  (4»d+tt/2)  f  (<j>D)  .  The 
geometry  arguments  also  apply  for  the  magnetic  field  as  well 
leading  to  the  statement 

HaB2(pD'Z+Z'  '  *d+7t/2)  "  HaB2(pD,Z+Z’'*D)  aaP'^'z  (61) 

We  now  collect  the  appropriate  results  presented  by  Banos, 
change  to  the  "D"  notation,  and  use  Equations  60  and  61  to 
obtain 
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H 


pB2 


■  „  /SB  1  3W22^ 


H$B2  =  S  sin(J)D  l  3z 


a2w 
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3Pt 


0  3  B 

HzB2  ■  S  coS*D  w 


3  o 

EpB2  =  Q  3?3p: 


(°22  -  C) 


E«B2  "  HPB2  *  HzB2  0 


’2B2 '  °(s7 +  )  (Gl2~  c) 


where 


H^B2  “  "S  Tp  (G22  "  C) 


4irk? 


(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 
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A  =  G22  “  G21  +  k2V22 


B  =  G22  "  Si  +  U22 


C  =  G21  “  22 


-1  ik2R 
G22  *  R  e 


G2l“Ve 


R  -  (4  *  I*’*')2)' 

R+  =  (p2  +  (z+z')2)1' 


U-n  »/ 


“Y2 (z+z' ) 


22  o  V*2 


«  -Y,Cz+z') 


V22  *  f 


•  k:Y,+krY, 

O  2  112 


JQ(XpD)XdX 


JQ(XpD)XdX 


.-■7 


<*>  -Y2(z+z') 

f  2(Y.-Y,)e  2 

j  - 5 - J  (Xp  )XdX 

°  ^Y1+kjy2  0  ■> 


29 


2  2  1/2 
r,  ’  U2-k2) 

(91) 

The  variety  of  branch  cuts  that  have  been  used  to  define  Y^ 
are  discussed  by  Banos.  We  are  now  in  a  position  to  utilize 
Equations  50  and  51  to  obtain  E^Cr^r'  )  and  H2 (£/£*)  and  then 
in  turn  to  use  Equations  39  and  47  to  obtain  the  desired 
quantities  g2 (r,r' )  and  g2(r,r').  Following  this  procedure, 
we  obtain 

(92) 

K2  (r ,r ' )  =  g0  (r , r 1 )  +  (r,r*  ) 

(93) 

where 

4irk2g 

0  (r,r ')  *  g8  a.  +  Gft(fl(fa.  a.  +  g„  a  a 

8 -  »PP  PD  PD  P'fxp  4>d  4>d  Spz  pd  z 

+  G8zpVpd  +  G8zzVz  0=o's 

(94) 

4’*3  <£'*’>  »  k8p*V*D  +  ***>\\ 

+  Kfizi*zhn  +  KW*zVz  6*°'s 

U  L 

(95) 

32g22  2 

Gopp  3pD5  +  k2G22 

(96) 
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(97) 


G._  =  — 


o4mJi 


1  ^73  2 

P~  T22~  +  k2G22 


D  3pr 


32G 


22 


aopz  3pD3z  ozp 


3  .2, 


22  +  klG, 


-o-  3z2  ‘2  22 


spp 


3p 


S‘  (-°21+k2V22)  +  k2(-G21+U22) 


2  21  2 


G.,,  +U„~ 
21  22 


s4>4> 


G3P2  *  Jp^5z  (-G21+klV22)  *  -G 


szp 


szz 


(l7  +  k2)  ("G21+klV22) 


-3G 


22 


-K 


o  p  3z  o4»p 


3G 


22 


-K 


oz<$  9p  o$z 


aw 


Ksp<J>  *  3z(G2l“U22)+  Pp  3p 


22. 

D 
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(98) 


(99) 


(100) 


(101) 


(102) 


(103) 


(104) 


(105) 


(106) 


(107) 


Ks$p  =  3z  ('G21+U22)"  3od2 

KS2*-  3^(-G21+U22)  ,108> 

Ks*2  *  3^  (G21-^22)  (109) 

After  making  the  described  substitutions,  it  was  necessary  to 
use  the  coordinate  relationships 

=  a  :cos<f>  +  a  sini)  (110) 

x  u  y  u 


and 


-&xsin<fiD  +  aycos4>D 


(111) 


in  order  to  obtain  Equation  92  through  Equation  109.  Even 

though  the  explicit  representations  for  a  and  a^  given  by 

PD  u 

Equations  110  and  111  were  used  to  obtain  the  expressions 
for  §2 (r,r*)  and  K2(r,r'),  it  may  be  more  convenient  to  use 
the  final  expressions  with  alternate  representations  for 
the  unit  vector  as  expressed  in  the  "0"  coordinate  system. 
From  Equation  56  it  follows  that 


(112) 


and  from  the  properties  of  the  orthonormal  set 
it  follows  that 


V 


t 
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X 


(113) 


a. 


p-p  • 


We  now  note  that  the  K  (r,r* )  and  S<,(r,r')  derived  in  this 
section  are  identical  to  the  corresponding  quantities  pre¬ 
sented  in  Reference  1  which  were  derived  by  a  different 
method.  To  complete  the  comparison  with  the  material  presented 
in  that  reference,  we  note  that  by  simply  performing  the  deriva¬ 
tives  indicated  in  the  representations  presented  in  that  ref¬ 
erence 


47rSo(r,r' )  =  (|+k"2w)  G22  (114) 

and 

4irK0(r,r')  «  VG22  x  I  (115) 

it  follows  that  these  expressions  are  identical  with  and 
£o  presented  in  this  section.  Finally,  we  note  that  the 
significance  of  the  representations  presented  in  Equations 
114  and  115  will  be  seen  in  a  subsequent  section  when  we  are 
concerned  with  the  limit  as  r  approaches  the  surface. 
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IV.  EXPLICIT  REPRESENTATIONS  FOR  E (r)  AND  Hf(r) 

FOR  THE  TWO-LAYERED  PROBLEM  OF  INTEREST  1 

In  this  section  we  will  present  the  same  results  that 
were  presented  in  Reference  1;  however,  we  will  rewrite  those, 
expressions  in  a  form  where  it  is  possible  to  readily  make  a 
very  important  observation.  We  will  restrict  our  attention 
to  the  case  where  the  incident  field  is  a  plane  wave  given  by 


Ei  *  eEoe 


ik2fio'£ 


1116) 


with  e  being  a  unit  vector  along  the  polarization  direction 
and  hQ  being  a  unit  vector  in  the  propagation  direction. 
Rewriting  the  expressions  for  ET  and  given  in  Reference  1< 
we  have 


*  -  f  /  ik.n  .r\ 

E0  &!•(£)  =  (6-hJh*  +  (§•$)  v<J>  +  2  H  *  °  ~j(fio.az)  (S2xh)  (117) 


v;1 


Ht(£)  -  -(e-h)j^v*  +  2  (fi0-  az>  (a^hnj  '118) 


where 


(119) 


ik,A  *r  ik.A  •£ 

*  -  •  2  °  "  +  RHe  2  oR 


(120) 


ik,A  *r  ik.A  *r 

$  -  e  2  °  “  +  Rye  2  oR 


(121) 
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ftrtT9  =  ft  -  2a  (S  •  ft  ) 
OR  O  2  Z  O 


(122) 


h  = 


A  A 

n^xa 
o  z 

IVM 


(123) 


v  =  h  x  n. 


(124) 


-k2(az.  ^-{kj-k^1 

hv  "o)  2] 

|l/2 

“k2(W+{*l“,t2 

'T'2 

(125) 


-kl(az*  fto)_k2^ 

2  2 1 
k  -k* 
1*1  *2 

i-(a  •  ft  )2' 

L  Z  O  J 

r 

fk5^kr 

r i  k2 

1-(S  •  ft  )7 
z  o 

yji 

(126) 


We  can  now  discuss  two  important  limiting  cases  of  E^,  and  H  . 
One  is  the  limit  as  the  conductivity  of  medium  1  becomes 
infinite  (k^-**”)  ,  and  the  other  is  the  limit  as  the  incident 
angle  becomes  grazing  (S^.ft^-’-o).  in  the  infinite  conductivity 
limit ,  we  have 


lim 

k^-*-« 


RH  =  -1 


(127) 


and  in  the  grazing  angle  limit,  we  have 


(128) 


j 
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-1 


(129) 


lim 

(S  -fiJ-o 
z  o 


*H  * 


lim 

(a  *n  )-*-0 
z  o 


*V 


=  -1 


(130) 


The  difference  in  sign  in  the  two  limits  expressed  by  Equations 
128  and  130  leads  to  considerably  different  behavior  for  the 
total  fields  in  these  two  limits.  An  important  case  to  con¬ 
sider  is  the  situation  where  the  observation  point  above  the 
ground  is  considerably  less  than  the  free-space  wavelength 
so  that  k2z«l.  For  this  case,  we  choose  a  convenient  coor¬ 
dinate  system  so  that  £  =  0  .  We  can  approximate  both 
ik2fi  *r  ik2fi  *r 

e  and  e  by  1,  so  that  under  this  restriction 

on  observation  height  we  have 

lim  ip  *  0  (131) 

k^-*-» 


lim  <t>  -  2  (132) 

k-j^00 


lim  <p  =»  0  (133) 

(S  -ft  )-0 
z  o 


lim 

<a  -ft  >• 


<P  =  0 


(134) 


Using  Equations  131  through  134,  we  have  under  the  restriction 
k2z<<l  the  following 
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(135) 


lim  e'1E_(0+) 
k1-»  o  T 


2{e.$)[0  +  (ft  •  a  )  (axh)  ] 
o  z  z 


lim  Z  E"1HT(0+)  =  2(£*v)h  (136) 

k^-«°  2 


which  after  algebraic  manipulation  can  be  written  as 


lim  e“1E(T1(0  +  )  =  2  (i  (137) 

k1'*-oo  °  1  2  Z 


lim  Z.E^H  (0+)  =  2  [  (ft  x£)*h]h  (138) 

k^oo  2  0  1  o 

Equations  137  and  138  are  readily  interpreted  as  representing 
the  expected  results  for  a  perfect  conductor,  i.e. ,  a  doubling 
of  both  the  normal  incident  electric  field  and  the  incident 
tangential  magnetic  field  with  all  other  components  vanishing. 
The  grazing  angle  limits  are  given  by 


lim  E’1E_(0+)  =  0 
° 

lim  z  E'1H_(0+)  =  0 
(h  •  S  )-*•“>  2 

O  Z 


(139) 


(140) 


The  implications  of  the  extremely  different  limits  of  the 
perfectly  conducting  limits  and  the  grazing  angle  limits  are 
important  to  note.  For  the  purposes  of  the  numerical  calcu¬ 
lations  performed  in  this  report,  we  note  that  the  induced 
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current  density  on  an  object  above  a  lossy  half-space  or  above 
a  perfectly  conducting  ground  will  be  significantly  different 
for  grazing  incidence  excitation.  More  important,  questions 
having  significant  implications  arise  as  a  result  of  the 
difference  between  these  limits.  An  immediate  question  is 
to  what  extent  will  a  multilayered  model  of  the  earth  and  its 
atmosphere  reduce  the  ground  coverage  of  an  EMP.  Another 
question  is  how  is  this  effect  accounted  for  in  one  of  the 
dominant  EMP  coupling  mechanisms,  i.e.,  the  near  grazing 
excitation  of  long  lines  above  the  earth.  Finally,  how  do 
we  interpret  past  ground-based  system  level  tests  where  the 
intended  exposure  was  grazing  incidence? 


V.  PRESENTATION  OF  THE  MAGNETIC  FIELD  INTEGRAL  EQUATION 
(MFIE)  AND  THE  ELECTRIC  FIELD  I NTEGRO- D I FFERENT I AL 
EQUATION  (EFIDE)  FOR  THE  STRATIFIED  MEDIUM  SCATTERING 
PROBLEM 


The  essential  details  in  arriving  at  the  MFIE  and  the 
EFIDE  for  the  stratified  medium  scattering  problem  are 
related  to  the  examination  of  the  limit  as  the  observation 
point  variable  r  approaches  the  surface  of  the  embedded 
scatterer.  Specifically,  we  are  interested  in  the 


limit  n  (r) 
r  -*-r  o  — 
—v  — 


x  Sn 


x  H„  (r  )  where  the 

— — V 


meaning  of  nQ(r) ,£,  and  ry  will  be  explained.  The  expressions 

for  <£v>  and  Hj,  (rv)  are 
s  s 


£Ns(£v)  =  ETN^(rv)  +  /gNs(rv,r')- J(r')dS'  (141) 

and 

^g-v*  =  — TNg  (£V)  {142> 

Equation  141  is  Equation  37  with  the  r  variable  relabelled 

r  ,  and  Equation  142  is  Equation  46  for  the  case  l  -  N  ,  again 
*  s 

with  the  r  variable  relabelled  rv»  The  reason  for  relabelling 
the  observation  variable  is  to  emphasize  that  rv  is  allowed 
to  vary  over  a  three-dimensional  volume  as  opposed  to  a  two- 
dimensional  surface.  We  reserve  the  variable  r  to  represent 
a  variable  related  to  ry  but  which  is  restricted  to  the  sur¬ 
face  of  the  scatterer.  In  general,  this  relationship  is 
expressed  as 


Ev  58  £  +  a 


(143) 
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and  we  are  interested  in  the  behavior  of  Equations  141  and 
142  as  |g[  =  q  approaches  zero.  For  many  situations  of  interest 
we  can  argue  that  the  angle  at  which  rv  approaches  r  does  not 
change  the  resulting  equations  and  conveniently  we  choose 


■  a  =  (144) 

where  ftQ(r)  is  the  outward  normal  to  the  scattering  surface 
at  point  r.  For  cases  wnere  Equation  144  is  not  applicable 
(e.g.,  r  corresponds  to  a  tip  or  edge),  it  is  still  true  that 
the  stratified  medium  problem  needs  to  be  treated  in  exactly 
the  same  manner  as  the  standard  free -space  MFIE  and  EFIDE 
limits.  Whether  or  not  we  can  employ  Equation  144,  we  will 
assume  that  a  unique  relationship  can  be  specified  between  any 
point  r  on  the  surface  and  a  point  £v  in  the  proximity  of 
the  surface  which  is  described  by  Equation  143  so  that 


lim  r  =  r 

lah°  ~v 


We  can  now  express  the  limits  of  interest  as 


li”  ao<£>  X  -  Urn  0o(r)  x  ETOjs<£v> 

+  lim  n  (r)  x  /  G„  (r  ,r ' ) • J(r ' )dS'  (145) 
-  o  -  J  =N„  -v  - - 


S  s 


lim  ft  (r)  x  H  (r  )  »  lim  ft  (r)  x  Hm„  (r) 
q-0  °  “  “N  -v  o  -  ~TN„  - 


■Ng  -v  ”6  -o 


+  lim  ft  (r)  x  f  K  (r  ,r')*  J(r')dS*  (146) 
q^O  °  ^  _„s  -v  - - 
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Using  the  continuity  of  the  fields  in  the  region  exterior  to 
the  scatterer,  the  boundary  condition 

n  (r)  x  E  (r)  =  0  (147) 

O  —  —  N  — 

S 

the  previously  cited  definition, 

ft_ (r)  x  H  (r)  =  J(r)  (148) 

°  _  -Ns  -  - 

and  -the  new  definition 

Ve)  x  Stn  <£>  =  Jt(£)  (149) 

s 

we  can  write  Equations  145  and  146  as 


-fi0(r)  x  Etn  (r)  =  IEL 
s 

J(£)  *  J^r)  +  iHL 

where 


(150) 

(151) 


V-’  xf  s« 


J(r '  )dS 


(152) 


ft^(r)  x  f 


(r  ,r') *J(r')dS' 

"V  —  — 


(153) 


The  limits  in  Equations  152  and  153  will  be  examined  by 

decomposing  GN  (£v»£')  and  gN  (£v>£')  into  parts  that  contain 
s  s 

a  singularity  and  the  remaining  nonsingular  parts. 
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The  singular  parts  will  be  expressed  in  such  a  manner  that 

the  limits  appearing  in  Equations  152  and  153  will  be  seen 

as  being  previously  evaluated  for  the  derivation  of  the  EFIDE 

and  the  MFIE  in  a  vacuum  environment.  To  accomplish  the 

desired  decomposition/  we  write  the  solutions  to  Equations  40 

and  41  for  the  case  SL  =  N  as 

=s 

li  (r  .r' )  =  E^ (r  » r ' )  +  E*(r  ,r')  (154) 

s 


and 


r') 


(155) 


The  quantities  E^(rv,r' )  and  H^(rv,r' )  not  only  satisfy  Equa¬ 
tions  40  and  41  for  £vevN  •  but  they  also  satisfy  these  Equa- 

s 

tions  considering  that  VN  constituted  all  of  space.  They  con 

s 

sequently  also  satisfy  the  radiation  condition  at  infinity. 

The  subscript  notion  U  is  introduced  to  indicate  that  all  of 

space  uniformly  consists  of  material  having  parameters 

UN  and  eN  .  With  this  understanding  of  the  meaning  of  E^ 
s  s 

and  H^,  it  is  possible  to  present  explicit  expressions  for 
these  quantities.  Instead  we  present  explicit  expressions  for 

3 

s0'v£'>  *  T,  <156) 

i*l 


and 


Su 


<£v'£'> 


3 


*  E  »i(v£,,ai 

i*l 


(157) 


which  because  of  the  described  meaning  of  E.  and  H.  can  be 

— 0  — u 

rewritten  as 

(iwyN  \ 

-^)(l^Vv)\(VRv) 


where 


Bu  <£*'£')  =  (41I) 


3  ,  /k„  ,  R  \  =  R~ 

M  Ns  v)  v 


=  R-^V* 


R  =  I  r  -r 'I 
v  I  — v  —  I 


Combining  Equations  39,  47,  and  154  through  157,  we  have 


<£*-£•>  *  §u(W  +  We'> 

s 


gN  <EV'E'>  *  Bu<Ev'E')  +  5r<Ev,E’) 

s 


with 


W'£’>  -  E5R'W>ai 

i«l 


SrIEv-E')  -  E«Jlrt,r'lli 


where  sufficient  background  has  been  presented  to  allow  the 
meaning  of  Ei(r,.,r')  and  H*(r_.,r’)  and  in  turn  G„  and  K„  to 


mmmm 


receive  further  elaboration.  ER  and  HR  are  solutions  to 
Equations  40  and  41  for  £  «  N_  with  the  6  function  source 
eliminated.  The  fields  E^  and  H*  satisfy  the  inhomogeneous 

—"U  —R 

equations,  and  the  purpose  of  ER  and  HR  is  to  cause  the 

fields  E^  and  given  by  Equations  154  and  155  to  satisfy 
s  s 

the  boundary  conditions  described  in  Appendix  A.  From  this 
understanding  of  the  nature  of  the  QR  and  KR,  we  can  conclude 
that  there  are  no  singularities  to  be  considered  in  the  limit 
as  q  approaches  zero  in  terms  containing  these  dyadics.  This 
is  expressed  explicitly  by  combining  Equations  152,  153, 

161,  and  162  to  obtain 


I  =  I  +  I 
-EL  iEU  -ER 


(165) 


and 


where 


T  as  T  +  X 

— HL  -HU  ^HR 


isR  *  x  f  §r{£'£.'  )  •£(£') ds 

S 

£)  x  / £R(r,r').J(r' 


“  fto ~  »R 


)  dS  * 


(166) 


IeU  =  lim  X  f =U(-v'-’ )  *-(£’ )dS'  (167) 


IHU  =  lim  nQ(£)  x  J Ky (£v,£' ) • J(r' )dS'  (168) 


(169) 


(170) 
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The  essential  behavior  of  the  limits  in  Equations  167  and  168 

2 

is  independent  of  whether  kN  is  real  or  complex.  For  the 
2  s 

case  where  k^  is  real,  these  limits  have  been  given  consid- 
s 

erable  attention  for  the  derivation  of  the  free-field  EFIDE 
and  MFIE.  Using  the  standard  results  for  these  limits,  we 
have 

k  *  ,ds; 


Ieu  *  ao(£)x 


£v^ 

(171a) 


or 


Ieu  =  fto(^)x 


(l  *  -  )dS'] 


(171b) 


and 


Iho  *  fto< 


n 


dS' 


+  *  *<£> 


^v=I 

(172a) 


or 


—HU 


fcQ(r)  x  f  (4Tr)_1^GN  ^kN  'K'j  X  J(r')jdS' 


x  4? 


(172b) 
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where  ft  is  the  solid  angle  subtended  by  the  surface  at  r. 

The  usual  expressions  for  the  limits  are  given  in  Equations 
171b  and  172b;  however,  the  real  meaning  is  better  understood 
by  referring  to  expressions  171a  and  172a.  It  should  be  noted 
that  none  of  the  integrals  in  Equations  171  and  172  are 
principal-value  integrals.  Finally,  we  refer  to  Reference  8 
for  alternate  expressions  for  Equation  171. 

Combining  Equations  150,  161,  165,  169,  and  171,  we  have 
for  the  EFIDE 

-ft0(£)  x  E^fr)  =  fi0(r)  x  |(j  ♦  k5^‘)j(lr)s jV*)--  >dS' 

♦  «„<£>  *  /[*.<£.£•>  - 

(173) 

Combining  Equations  151,  162,  166,  170,  and  172,  we  have  for 
the  MFIE 

^  ‘  =  Jt(£)  +  ftQ(£)  x  (4tt)-1^7Gn  ^kN  ,R^jx  J(r')dS' 

+  hQ(r)  x  (r,r  ’ )  -  K^r,^)!  -jfr'JdS' 

(174) 

There  are  two  important  points  to  note  about  the  representations 
given  by  these  equations.  The  first  is  that  this  decomposition 
is  not  just  formal  as  can  be  seen  by  observing  Equations  92, 

93,  114,  and  115  for  the  two-layered  medium.  The  second  is 
that  the  numerical  treatment  of  these  equations  can  be  facil¬ 
itated  by  using  these  representations.  For  the  MFIE  numerical 
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solution  of  the  two-layered  scattering  problem,  we  use  this 
decomposition  and  pay  particular  attention  to  the  treatment 
of  the  integrable  singularity  that  exists  in  the  first  integral 
appearing  in  Equation  174. 


VI.  EXPLICIT  SYMMETRY  ARGUMENTS  FOR  THE  TWO  LAYERED  MFIE 


The  actual  surface  that  was  numerically  treated  had  a 
plane  of  symmetry  perpendicular  to  the  ground  plane,  and  this 
allowed  us  to  reduce  the  computer  time  to  solve  the  problem 
numerically.  For  the  purposes  of  illustration,  we  will  choose 
our  coordinate  system  so  that  the  symmetry  plane  of  the  body 
corresponds  to  the  plane  x  =  0;  however,  the  final  results  will 
appear  in  a  form  that  is  invariant  to  this  coordinate  choice. 
First  we  write  Equation  174  as  specialized  to  the  two-layered 
problem  discussed  in  Section  III  as 

f  (ft) J(r)  =  J^r)  +  /  S(r,r').  J(r')dS'  (175) 

S 

where 


f<£!)  -  1  -  £ 


(176) 


and 


Q(r,r')  -  ftQ(r)  x  g2(r,r' )  (177) 

and  alternate  explicit  representations  of  £2  (r,r')  are  obtained 
by  combining  Equations  93,  95,  106  through  109  and  either  104 
and  105  or  115  depending  on  our  choice  of  representation  for  the 
KQ(r,r')  portion  of  £2(r,r' ).  For  the  choice  of  Equation  115 
to  represent  Ko(r,r’ ),  all  necessary  symmetry'  arguments  have 
been  presented  in  Reference  1,  and  we  could  make  use  of  those 
arguments.  Instead,  it  provides  little  more  effort  to  con¬ 
sider  Equations  104  and  105  for  the  representation  and  thus 
presents  a  self-contained  derivation. 
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We  now  write  Equation  175  as 


f(fl)J(r^)  = 


+ 


g(r+ 


,  r ' ) .J(r’)dS' 


r ' )  . J(r ' ) dS ' 


(178) 


and 


f(G)J(r  )  =  £T(r_) 


+ 


f B(r  ,r|) . J(r|)dS' 
s+" 


+ 


Jd(t  ,r ' )  . J(r')dS‘ 
S_  “  ” 


(179) 


where  r+  and  i'+  describe  points  having  positive  x-  coordinates 
while  r  and  r'  describe  points  having  negative  x- coordinates. 
The  integration  over  S+  is  over  the  surface  having  positive  x 
coordinates,  and  the  integration  over  S_  is  over  the  surface 
having  negative  x- coordinates.  The  explicit  relationship 
between  symmetric  points  is 

r  =  R  .  r  (180a) 

~x  + 


where 


r'  *  R  .  r'  (180b) 

— -  =  x  + 


B  *  I  -  (181) 

~x  x  x 
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We  now  rewrite  Equations  178  and  179  after  operating  on 
sides  of  Equation  179  with  Rjj  to  obtain 


f(fl)J  (r+) 


=  Jj(r+)  +  /  D++(r+,r;) .j+(r;)dS* 


+  f  s+- •J“(?;)ds’ 


(r+) 


f  +  ■ 

=  Jr (£+)+  {  U-+ (£+'£|.> •£  (?+)dS' 


+  f  B__Cr+,r|) .j"(r;)dS' 


where 


J+(r+)  -  J(rf) 

j+(£+>  ■  £(r;) 

=  V-V 


J  (r.)  =  R  • J(B  *r  ) 
—  x  —  -x  -+■ 

£‘<n>  -  sx-i(Sx-r;> 


both  x 

(182) 

(183) 

(184a) 

(184b) 

(185) 

(186a) 

(186b) 
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^(r+)  -  Sx*  } 


(187) 


B ++(r+,r;)  =  B(r+/r;>  (188) 

B+_  <?+'£+)  =  B(r+,Bx*r|)*Rx  (189) 

8_+(r+,r;)=  gx- D(gx-r+,r;)  (190) 


-  Bx-g<gx-r+.gx-£i>-Sx  (191) 

and  in  order  to  obtain  Equations  182  and  183  we  substituted 
the  unit  dyadic,  I,  expressed  as 

I-Bx'8x  (192) 

to  obtain  the  relationships 


8<F+'El>-  i(£l)  =  B(r+,r^).gx  .  gx.  J(rM 


(193) 


Q(r_,r' )•  J(r^)  =  (£_'£.' )  •  Sx  *  Bx‘ £(£ ' ) 


(194) 


which  in  turn  are  utilized  in  those  equations.  If  the  following 
relationships  are  true. 
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§++<*+' 2: +>  =  D__  (r+,r|)  =  &(r+,r|) 


(195 


and 


D+_(r+,r|)  =  B_+ (£+*£!>  =  §<£.+  /£+> 


(196 


then  by  first  adding  and  then  subtracting  Equations  182  and 
183/  we  obtain 


ff(B) 


V-V  *  irs<f+>  +  /  Bs^'U'-is'U*35'  <197 


and 


f«2) 


-  irA(£+>  +  /  JA(?;)ds' 


(198 


where 


i-S  ^-f}  =  f 


J  (rf)  +  J~ (r+) 


(199a 


-s^^V  ”  ? 


;j+(r;)  +  J'(r^) 


(199b 


irs'-r;'  '  I  2j<E+l  +  £<£+: 


(200 


5Al£+)  -  5  J+(£+)  -  j"<r+) 


(201 


a*(n>  - 1  j"(£;>  -  £'(r;> 


(202; 
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(203) 


W-V  =  t 


J^(£+)  -  :V£+> 


2s(?+,r_J_)  =  A(r+,r|)  +  §(?+,r_J.) 


(204) 


2a (£+,£+>  “  ^(?+'£|}  “  +  (205) 

and  the  advantage  of  Equations  197  and  198  over  175  will  be 
explained.  Following  that  explanation,  a  proof  that  Equations 
195  and  196  are  true  will  conclude  this  section. 

The  essential  property  that  the  set  of  Equations  197  and 
198  have  is  that  only  half  of  the  symmetric  surface  requires 
numerical  treatment.  The  resulting  matrix  that  gets  generated 
for  each  of  the  Equations  197  and  198  is  half  as  large  as  the 
matrix  that  would  get  generated  by  the  straightforward  numerical 
treatment  of  Equation  175.  Because  matrix  inversion  time  is 
proportional  to  the  cube  of  the  matrix  size,  there  is  a  factor 
of  four  savings  in  inverting  two  matrices  that  are  half  the 
size  of  a  larger  matrix  since 

2(!)3  -  T  »3 

Regarding  the  matrix  generation  time,  there  is  a  forced  savings 
caused  by  seeking  the  relationships  expressed  by  Equations  195 
and  196.  The  reason  there  is  a  forced  savings  is  that  the 
straightforward  treatment  would  effectively  require  that  each 
matrix  element  be  viewed  as  receiving  a  contribution  from 
g++,D _ ,D+_,  and  D_+,  while  Equations  195  and  196  show  that  there 

are  only  two  independent  contributors.  Returning  to  Equations 
197  and  198,  we  numerically  solve  for  Jg  and  JA  which  in  turn 
yield  the  desired  physical  quantities  J(r+)  and  J(r_)  by 
employing  the  relationships 
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L 


H 


£<r+)  =  Js(£+)  +  JA(r+) 


(206) 


J(r_) 


R  • 
=x 


Js(r+)-3A(r+) 


(207) 


Returning  to  the  proof  of  the  validity  of  Equations  195 
and  196,  we  identify  the  explicit  meaning  of  these  quantities 
as 


S++ <£+'£;>  - 

h  (r.)  x  K (r  ,r ' ) 

O  —  "— ■+■  —  f 

(208) 

!..(r+,r;)  =  Rx- 

nQ(r_)  x  K(r+,r^)|  *  §x 

(209) 

B+-<£+'£i>  3 

iQ(r+)  x  E(r+,r^)|  .  §x 

(210) 

g_+<I+»r|)  *  §x 

•  [(^o (£_)  x  £(r_,r;)| 

(211) 

and  the  definition  of  K(r,r')  is  given  by  combining  Equations 
93,  95,  and  104  through  109.  For  the  purposes  of  this  section, 
we  write  the  essential  form 


K(r,r')  =Kp^(r,r')a  (r,r')S  (r,r’) 

+  Knrfi(r,r')a.  (r,rf)a  (r,r’) 

P<f> - ?D - PD 

+  v^’^zV^') 

+  (£,£')  (£,£')  az  (212) 
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with  the.  essential  property  extracted  from  those  equations 
that 


Kag 


Kas 


(213) 


and 


Kag<£f'£l>  B  Kag(£-'I+>  (214) 

with  a  and  8  taking  on  the  notation  p,$f  and  z  to  generate  the 
four  terms  appearing  in  Equation  212.  Equations  213  and  214 
follow  by  noting  that  according  to  Equations  105  through  109 

Kag(£'£'>  =  K0a8  (l£-r'|)  +  Kgag  ( |  p-p  *|  ,z+z’)  (215) 


The  more  difficult  aspect  of  proving  Equations  195  and  196  is 
related  to  the  dyadic  character  of  these  equations.  To  complete 
this  aspect  of  the  proof,  we  will  introduce  the  following  dyadic 
identities  which  will  later  be  proved. 


fix  '  (s«)"  -(bx-*)*  (fix's) 

(216) 

fix  •(*«)-  -(fix'S)x(fix'£) 

(217) 

(^s)'fix  -  »x(s-fix) 

(218) 

Using  these  equations  as  well  as  Equations  212, 
we  can  write  Equations  208  through  211  as 

213,  and  214, 

B++<r+.s;»  -  Vr+>  x  jKp«<r-i.-Ei> spd(j+'I+>s4d 


-  %p(r*.r;)»  (r+,r;>a  (r+,r;> 


♦  K-lr+.r;)*^  (r+,r;> 


♦  k*,  <*+.*;>  **D<r+.ryaz| 


(219) 


B~  <£+»*;>  -  -Bx*A0(r  Jx  Kp<j) <r+,r;)  -gx 


+  [§*v£-'£-)]  [h^.^% 


♦  *,♦<£♦•£;>  (fix-*,)  [SD IE..*:) -b*; 


+  %zi£+-r;>[Bx-  **D(£-'£-,lt4*‘  Bx>j 


(220) 


8+- <£♦•£'♦>  -  »0(£+)  x  jKpp(£t.»LI»0_(r+rrll 


*♦„<£+•£:.>  -8x 


*  K*p(r+.r:)»  <r+.r •_) 


sPD(r+'El)-8x 


+  *M  Ir+.r-J  Sz 


V <£»■£!> -Sx 


*  %g(£+.r;)*p -g )  (221) 
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B.+ <r+-ri»  * 


'  Bx-flo<£.>]  x|K0«(-r+'£:)[Sx’ipD(£.'-ril 


)  a.  (r  ,r* ) 

4>d  ~ 


+  %p<r+.£:)[?x-^o<r_,r;)]  SPd<£_.£;> 


+  K2t  '£»'£:>  <Bx'iz>\<£_'Zl> 

*  Kn  «£+'£:>  [ftc-SD«£.'-r4>]4,| 


(222) 


The  equations,  which  when  substituted  into  Equations  219 
through  222  lead  to  the  desired  results,  will  be  presented. 
First  we  present  the  following  representations  for  where 
x  and  x’  are  positive  D 


a  (r.,r!)  * 
PD  -+  -+ 


(x-x* )£ +  (y-y’ )a 


-(x-x’>  a  +  {y-y ' )a 
a  (r  ,r’_)  =  - -s - i 


D  - 
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(x+x')a+(y-y’)S, 

X  > 


-(x+x')S  +(y-y’)S 

K  «£-<£;> - r - 1 


P_  =  {x-x')2+(y-y’)2] 


p+  *  (x+x' ) 2 


+  (y-y')2] 


(223) 


(224) 


(225) 


(226) 


(227) 


(228) 
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By  simply  applying  the  indicated  operations  on  the  preceding 
equations,  it  follows  that 


5x-ap  (£  ,r ' )  =  4  (£_,£')•£  «  4  (r  ,r;) 

x  PD  PD  PD 


(229) 


Vr')-R  *  4  (r  ,x\) 
md  mD  “ 


(230) 


Bv*  a  (r  ,r ' )  »  1  (r.,r') 

-x  pD  -+  PD  -+ 


(231) 


Next  using  the  definition 


a.  (r,r ' )  =  4.  x  an  (r,r*  ) 
<frD -  2  PD - 


as  well  as  Equation  216  and  the  fact  that 


R  -a  =  a  *R»  4 
=x  z  z  -x  z 


(232) 


we  find  that 


B  *a.  (r  ,r  ')  =  a.  (r  ,r')*B  «-a.  (r.,r!)  (233) 

-x  <PD  --  4>d -  **  ♦D-+-  + 


(r+,£')-Rv  -  “4  (r  ,r') 
♦D  +  “  x  *d  + 


(234) 


Bx*%D(r.r;) 


-4*  (E+'Sl5 


(235) 
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Hie  final  relationship  needed  follows  from  the  basic  symmetry 
of  the  body,  and  it  is 

gx.fio(rJ  =  nQ(r+)  (236) 

Substituting  Equations  229,  230,  231,  as  well  as  232  through 
236  into  Equations  219  through  222  yields  the  desired  results, 
i.e..  Equations  195  and  196. 

This  section  is  concluded  by  presenting  the  derivation  of 
Equations  216,  217,  and  218  as  promised  prior  to  the  use  of 
these  equations.  First,  we  decompose  general  vectors,  A  and  B, 
as  follows 


A  =  Ay  +  Axax 


(237) 


and 

S  -  Sj.  +  Bxax  (238) 

where  the  subscript  T  indicates  a  direction  orthogonal  to 

a  •  We  now  form  the  cross  product 
x 

A  x  B  =  A_  x  B_  +  A  a  x  B_  +  A_  x  B  a  (239) 

—  —  — T  — T  xx— T^T  xx 

and  note  that  the  first  term  on  the  right-hand  side  of  Equation 
239  is  directed  strictly  along  the  ax  direction  while  the  second 
and  third  terms  have  no  ax  directed  components.  Using  this 
observation,  it  follows  that 

Ex'  <i*5>  -  -VS-T  +  AxV  WBx*x 


(240) 


Operating  on  the  expressions  given  in  Equations  237  and  238 
and  then  forming  the  cross  product,  it  follows  that 


(Sx'i>  x  *  WVx^T-Wx  (241» 

Multiplying  this  equation  by  minus  one  and  substituting  the 
result  into  Equation  240  yields  the  desired  result,  i.e.. 
Equation  216.  Next  we  express  a  general  dyadic  Q  as 


S  *  £  Mi 

i=l 


(242) 


so  that  we  can  write 


3 

Z 

i=l 


A  X  g  =  J]  (AxEi> F± 


(243) 


Operating  on  this  equation  to  the  left  with  gx  yields 

3 


Ix* (Ax£> 


-  £ 

i=l 


Bx*  (Ax^) 


F. 
— i 


(244) 


any  by  utilizing  Equation  216,  we  obtain 

gx' (Axg)  *  Yj  ■(§x’-)  x  ( =x '  — i } - i 
i*l 


-(§x^  x  §x‘  Z 

i»l 


(245) 
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which  becomes  Equation  217  after  resubstituting  Equation  242. 
Writing 


(Axg)  -gx 


R 


=x 


b.  *  E 

i»l 


i*(t -  M 


(246) 


completes  the  derivation  of  Equations  216,  217,  and  218. 
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VII.  CODE  TESTING 


1.  INTRODUCTION 

The  incorporation  of  a  lossy  half  space  into  our  MFIE  code 
(NEC-2A) ,  described  in  Reference  1,  was  complicated  by  the  absence 
of  data  related  to  the  model  depicted  in  Figure  4  that  could  be 
used  to  quantify  the  accuracy  of  the  added  capability.  Special 
tests  had  to  be  devised  to  detect  the  presence  of  coding  errors. 

In  this  section  we  describe  the  steps  taken  to  establish 
confidence  in  the  added  code.  We  first  summarize  these  steps 
and  then,  for  the  interested  reader,  discuss  the  detailed  con¬ 
siderations.  At  the  end  of  this  section  we  discuss  some  of  the 
accuracy  considerations  for  the  numerical  procedures. 

•  As  can  be  seen  from  Equations  62  through  73,  the 
computer  code  would  require  the  derivatives  of 
U22'  V22  an<*  W22‘  Although  we  were  unable  to 
obtain  tables  of  necessary  quantities,  we  did  find 
tables  for  U22  and  V22  (Equations  89  and  90).  We 
therefore  developed  computer  codes  for  calculating 
these  integrals  and  validated  these  codes  over  the 
range  of  the  tables.  We  then  developed  the  subroutines 
that  would  be  used  to  calculate  the  U22,  V22  and  W22 
derivative  terms.  The  U22  and  V22  derivatives  were 
verified  by  direct  comparison  to  the  numerical  deri¬ 
vatives  of  the  terms  previously  calculated.  The  W22 
derivatives  were  validated  indirectly  through  relation¬ 
ships  between  U22,  V22  and  W22> 

•  New  coding  was  required  to  incorporate  the  new  terms 
into  the  matrix  of  influence  coefficients.  By 
examining  the  limit  of  a  perfectly  conducting  half¬ 
space,  we  were  able  to  verify  that  the  new  code 
properly  communicated  with  the  other  modules  of  the 
computer  code. 


62 


+(sin$psin<i0cosa0  +  cos4pcos«0)ey 
-s1n4pSine0ez 


Figure  4.  Description  of  the  incident  plane  electromagnetic 

wave.  Axis  x"  is  the  intersection  of  the  z,k  plane 
and  a  plane  perpendicular  to  k  at  the  origin. 
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•  In  the  limit  of  the  half-space  having  the  dielectric 
properties  of  f ree-space ,  the  code  should  produce  its 
free-space  answers.  In  this  manner  we  were  able  to 
validate  the  relative  sign  of  the  additional  terms. 

•  The  Lawrence  Livermore  Laboratory  Nuclear  Electromag¬ 
netics  Code  (NEC- (1, 2)  A.)  can  calculate  the  bulk  current 
for  a  special  case  of  the  situation  considered  in  this 
report.  By  integrating  the  longitudinal  current  densi¬ 
ties  of  our  own  predictions,  we  showed  that  our  results 
agreed  to  within  five  percent  of  those  produced  by  the 
NEC- (1,2) A  code  for  our  test  case. 

•  Our  final  test  was  a  comparison  to  experimental  data 
obtained  at  the  University  of  Michigan  for  a  situation 
that  could  presumably  be  modelled  by  our  code.  We 
made  comparisons  for  a  cylinder  in  free  space,  a 
cylinder  above  a  perfectly  conducting  ground  plane, 
and  a  cylinder  above  a  planar  lossy  surface.  Quali¬ 
tatively,  there  is  agreement  for  all  three  situations; 
quantitatively,  there  are  explainable  discrepancies 

at  resonances  for  the  perfectly  conducting  cases,  but 
the  discrepancies  for  the  cases  involving  a  planar 
lossy  surface  are  larger  than  expected.  It  is 
possible  that  the  quantitative  discrepancies  for 
the  lossy  planar  surface  are  due  to  differences 
between  the  numerical  model  and  the  physical 
situation  in  which  the  lossy  half-space  was  an 
inadequate  representation  of  the  lossy  material. 

We  will  now  amplify  the  above  summary  for  those  readers 
who  are  interested. 
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2.  VERIFICATION  OF  BASIC  QUANTITIES 


As  explained  in  the  beginning  of  this  section,  we  were 


au 


unable  to  obtain  tables  for  the  quantities 


22 

9P,  ' 


3U 


22 


3V 


22 


3z 


'D  "  3PD  ' 

3W22  3W22 

and  — =^-.  Instead,  we  developed  subroutines  to  compute 


3Pr 


3P 


U22  and  V22  and  verified  their  accuracy  through  comparison  to 
published  tables  (Reference  9).  Having  verified  the  accuracy 
of  these  subroutines  we  proceeded  to  calculate  the  partial  deri¬ 
vatives  of  U22  and  V22  by  applying  finite  difference  techniques 
to  the  results  provided  by  the  above-mentioned  subroutines. 

These  finite  difference  estimates  of  the  partial  derivatives 
were  in  turn  used  as  our  reference  for  validating  the  sub¬ 
routines  used  in  the  larger  computer  code. 

It  was  our  judgment  that  the  above  partial  derivatives 
could  be  more  economically  and  more  stably  calculated  if  they 
were  analytically  included  inside  the  integrals  defining  the 
required  terms.  We  therefore  coded  these  quantities  as  follows: 


„  .2  f 

5»d  J 


l ' vpn'  9 9 

1  D  k":Y2+k,Y1 


(249) 


3pD  V 


™22  -  -*»  ^*dXe~Y2  (z+z’}  j  (Xp  )  — ~— 

k  Y2+koY1 


(250) 


!^.-2  AVv->  LapD,  -CT^f^un, 


3P 


J i  (Xp_ ) \  (Y,“Y, ) X' 


D  /  k  Y2+k;Y1 


where  pp,  z,  and  z*  are  defined  as  in  Equations  54,  55  and  56 

k0  =  Wave  number  of  the  incident  wave  in  free  space, 
k  =  Complex  wave  number  of  the  wave  transmitted  through 
the  lossy  half- space;  k  has  a  positive  imaginary 
part. 

*^22  *^22  ^22 

The  accuracy  of  -r- — ,  — * —  and  —  were  checked  by 

9pd  '  az  apD 

comparison  to  the  previously  established  finite  difference 

aw 

results. 

"jd  ap 


a2w„ 


22  22 

and  - =—  could  not  be  checked  by  these  means, 

D 


We  indirectly  validated  the  quantities  by  showing  that  they 
satisfied 


5S2  „  1  8W22  _  .  8G21 

A  '  n“  a  _  “  “ 


3p 


a 


PD  3p, 


az 


+  k2f^2  +  8U22 


az 


3z 


(252) 
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3.  USAGE  OF  THE  LOSSY  HALF- SPACE  TERMS 


The  matrix  is  computed  as  the  sum  of  three  matrices: 
a  free-space  matrix,  a  matrix  due  to  the  presence  of  a  per¬ 
fectly  conducting  ground  plane,  and  a  correction  matrix  to 
account  for  the  discrepancy  between  a  perfectly  conducting 
half-space  and  a  lossy  half-space.  This  decomposition  is 
described  in  Reference  1. 

The  form  of  the  free-space  operator,  consistent  with  the 
actual  programming,  is 


n  (r) 
4  it 


(r-r') 


iklr-r' | 

x  i  ® - - (ik|r-r'  |-l)dA 

*  i _ •  |J  ~  ” 


(253) 


and  the  additional  operator  for  a  perfectly  conducting  ground 
plane  which  is  in  a  form  consistent  with  the  actual  programming 
is 


§*n 

4ir 


x 


(R;r-r‘) 
*2—  — 


x  I 


(ik  f R* r-r ' [ -1) dA' 


(254) 


where  by  R^a  is  meant  the  reflection  of  a  about  the  ground 
plane.  As  can  be  seen,  the  two  operators  have  the  same  form 
with  respect  to  the  variable  of  integration,  r' ,  thus  permit¬ 
ting  the  use  of  the  same  subroutines  for  computing  both  integrals. 

The  lossy  half-space  correction  term,  however,  cannot  be 
cast  in  this  form;  hence,  separate  subroutines  are  needed.  The 
lossy  half-space  operator  is 


n(r)  x  J 


/V  A 


a-  +  K  Aa  a. 
s$?  z  szq>  z 


A  A 


A  A 


+  Ksp<t>a4>DapD  +  Ksp4>aPDa<J>D 


dA  (255) 


where  K  .  ,  K  , ,  K  .  and  K  .  are  given  by  Equations  106-109, 

s$z  sz<p  sp<p  Stpp 

£  is  the  unit  vector  perpendicular  to  the  ground  plane, 
z 

£  is  the  direction  of  the  projection  onto  the  ground  plane 

PD  . . 

of  the  vector  'from 'r''*  to’ r  (if*  the  projection  has  a  zero 

length,  any  given  unit  vector  perpendicular  to  £  will  do) , 

z 


i  x  £ 

2  P. 


and  a . 

“  KD 

Fortunately,  the  perfectly  conducting  operator  (and 
therefore  the  free-space  operator)  can  be  cast  in  the  same 
form  as  the  lossy  half-space.  For  the  perfectly  conducting 
operator. 


K  =  K 
SZ$  Stpz 


-elk° “iz---' I (ik0 | R- r-r ’ |-1) pd 
_  "z -  D  (256) 


4  ir  |  R-r-r  ' 


and 


.^Jg-r-rV 


K  .  =  -K  . 

s4>  P  sp<}) 


(iko  IgjE’E'  |-l)(z+z' ) 


(257) 


4ir  t  R-  r-r' 

1  =z - 


By  defining  the  scalar  quantities  so  as  to  match  the  perfectly 
conducting  operator,  we  were  able  to  test  the  subroutines  that 
were  written  to  incorporate  the  lossy  half-space  terms.  These 
subroutines  produce  the  same  results  as  the  subroutines  that 
calculate  the  perfectly  conducting  operator. 

4.  FREE-SPACE  LIMIT 


As  discussed  in  the  preceding  subsection,  the  total 
operator  is  written  as  the  sum  of  three  operators.  For  the 


special  case  of  the  half-space  having  dielectric  properties 
similar  to  those  of  free-space,  the  total  operator  should  be 
similar  to  that  of  the  free-space  operator;  i.e.,  the  second 
•  *and  third-'operator-s  should  have  .similes  .magnitudes  but  opposite 
signs.  Having  verified  that  the  code  had  this  property,  we 
exhausted  the  available  checks  on  the  incorporation  of  the 
new  terms. 

5.  COMPARISON  TO  ANOTHER  CODE 

G.  J.  Burke  of  Lawrence  Livermore  Laboratory  provided 
us  with  predictions  by  the  NEC- (1,2) A  code  of  the  bulk 
current  induced  on  a  cylinder  in  the  vicinity  of  a  lossy 
half-space.  The  parameters  of  the  case  are  shown  in  Table 
1.  The  results  of  the  two  codes  are  shown  in  Table  2.  Our 
results  were  obtained  by  a  four-point  averaging  of  the  longi¬ 
tudinal  current  densities  at  a  given  distance  from  the  center 
of  the  cylinder.  The  discrepancy  between  the  two  predictions, 
less  than  5  percent,  can  be  attributed  to  the  large  zone 
sizes  used  by  both  codes,  although  certain  differences  in 
the  modelling  assumptions  could  also  contribute  to  this  small 
error. 

6.  COMPARISON  TO  EXPERIMENTAL  DATA 

Here  we  present  comparisons  of  our  code's  surface  current 
predictions  to  experimental  data  obtained  at  the  University 
of  Michigan.  The  first  three  figures  (Figures  5-7)  are  repro¬ 
duced  from  Reference  1;  they  are  included  as  indications  that 
the  experimental  equipment  is  measuring  the  quantity  being 
predicted  by  our  code.  The  remaining  three  figures  (Figures 
8-10)  show  how  the  experimental  data  for  a  cylinder  above  a 
lossy  planar  surface  compares  to  our  predictions  of  surface 
current  for  a  cylinder  above  a  half-space  having  the  same 
dielectric  properties  as  the  planar  surface.  Table  3  de¬ 
scribes  these  dielectric  properties.  For  each  graph,  the 
solid  curves  are  subjectively  smoothed  traces,  i.e.,  spikes 


TABLE  1.  CONFIGURATION  PARAMETERS 
FOR  COMPARISON  TO  THE 
LAWRENCE  LIVERMORE 
LABORATORY'S  NUCLEAR 
ELECTROMAGNETICS  CODE 
(L3NEC) 


Wavelength  = 

1  m 

Cylinder  radius  = 

0.025  m 

Cylinder  length  = 

0.5m 

Cylinder  height  above 

half-space  = 

0.1m 

j  Relative  dielectric  1 

constant  = 

8  +  iS 

TABLE  2.  PREDICTED  BULK  CURRENTS 


Distance  from 
midpoint  of 
cylinder  (m) 

MFIE 

l3nec 

Percent 

difference 

0.0 

0.00395 

0.00413 

4.6 

0.0714 

0.00370 

0.00380 

2.7 

0.1429 

0.00289 

0.00287 

0.7 

0.2143 

0.00151 

0.00144 

4.6 
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Figure  5.  Comparison  between  "Free  Cylinder"  experimental  results 

at  Unviersity  of  Michigan  (solid  curves)  and  MFIE  computer 


Figure  6.  Comparison  between  "Cylinder  Parallel  to  a  Perfectly  Conducting 
Ground"  experimental  results  at  University  of  Michigan  (solid 
curve)  and  MFIE  computer  code  (dots)  for  d  =  5a. 


inc 


Figure  7.  Comparison  between  "Cylinder  Parallel  to  Perfectly  Conducting 
Ground"  experimental  results  at  University  of  Michigan  (solid 
curve)  and  MFIE  computer  code  (dots)  for  d  =  1.5a. 
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TABLE  3.  COMPLEX  DIELECTRIC  CONSTANT  USED  IN 
THE  MFIE  COMPUTER  CODE  CALCULATIONS 
APPLICABLE  FOR  THE  UNIVERSITY  OF 
MICHIGAN  EXPERIMENTAL  DATA 


Actual  frequency 
of  experiment 
(GHz) 

kh  used  in 
the  code 

e2R/€o 

e2l/eo 

0.61 

0.7 

11.7 

15.6 

0.78 

0.9 

10.4 

13.2 

0-87 

1.0. 

9.6 

12.5 

0.95 

1.1 

9.0 

12.0 

1.00 

1.15 

8.7 

11.6 

1.04 

1.2 

8.6 

11.2 

1.09 

1.25 

8.5 

11.0 

1.13 

1.3 

8.1 

10.6 

1.17 

1.35 

7.9 

10.4 

1.22 

1.4 

7.8 

10.2 

1.30 

1.5 

7.3 

10.0 

1.48 

1.7 

6.8 

9.3 

1.74 

2.0 

6.3 

8.2 

2.08 

2.4 

5.7 

7.9 

2.43 

2.8 

4.6 

7.3 

2.78 

3.2 

4.3 

6.5 

3.13 

3.6 

4.2 

6.0 

3.47 

4.0 

4.0 

5.4 
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Figure  10.  Comparison  between  the  MFIE  computer  code  (dots)  and 
the  cylinder  parallel  to  a  finitely  conducting  ground 
experimental  results  at  University  of  Michigan  (solid 
curve)  for  d  =  1.5a. 


eliminated,  of  experimental  data  while  the  dots  are  the  values 
predicted  by  our  code;  see  Reference  1  for  the  original  experi¬ 
mental  curves. 


We  observe  that  except  at  resonance  the  agreement  between 
prediction  and  experiment  is  excellent  for  the  free-space  case 
and  the  case  of  a  cylinder  above  a  perfectly  conducting  ground 
plane.  The  fact  that  the  agreement  suffered  most  at  resonance 
is  not  surprising  since  this  region  is  most  sensitive  to  de¬ 
tailed  modelling  accuracy. 

The  same  cylinder,  source  and  probe  were  used  to  measure 
the  cylinder's  response  to  the  incident  wave  in  the  free- 
space,  perfectly  conducting  half-space,  and  lossy  half-space 
configurations.  Therefore,  we  will  use  the  level  of  agreement 
obtained  for  the  two  former  configurations  as  a  basis  for  com¬ 
paring  the  predictions  to  the  experimental  data  for  the  lossy 
material.  According  to  this  criterion,  although  there  is  gen¬ 
eral  qualitative  agreement,  the  quantitative  agreement  suffers 
by  comparison  to  the  earlier  efforts. 

This  disagreement  does  not  necessarily  invalidate  either 
the  experiment  or  the  computer  code.  An  alternative  is  the 
possibility  that  the  differences  between  the  experimental 
configuration  and  the  configuration  being  modelled  by  the  code 
are  significant.  Take,  for  example,  the  lossy  ground  plane 
used  in  the  experiment.  It  consisted  of  12  slabs  assembled 
into  three  layers  4  ft  x  4  ft  x  2*s  in.  supported  by  a  metallic 
ground  plane.  How  this  compares  to  a  half-space  is  in  question 
(Of  course,  if  we  believe  both  the  experimental  and  predictive 
results,  the  question  is  answered.) 

7.  ACCURACY  CONSIDERATIONS 

The  MFIE  is  a  direct  consequence  of  Maxwell's  equations 
for  the  problem  being  solved:  electromagnetic  scattering  from 
a  perfect  conductor  embedded  inside  a  layer  of  a  stratified 


medium.  However,  the  numerical  approximation  used  for  solving 
this  equation,  patch  zoning,  cannot  be  performed  with  complete 
precision  in  a  finite  time.  Besides  the  approximation  inherent 
in  the  technique  (piecewise  constant  solutions) ,  the  matrix 
elements  must  be  approximated  since  the  kernels  of  the  asso¬ 
ciated  inteqrals  cannot  be  integrated  in  closed  form;  in  fact, 
for  the  lossy  half-space  terms,  the  kernels  are  themselves 
integrals  which  cannot  be  integrated  in  closed  form. 

The  parameters  affecting  solution  accuracy  are: 

1.  Zoning  density 

2.  Quadrature  order  for  spatial  integration 

3.  Accuracy  of  the  terms  due  to  the  lossy  half-space 

The  first  item  has  been  programmed  as  input  variables; 
the  user  of  the  program  decides  their  values.  Our  choice  for 
the  second  and  third  parameters  was  dictated  by  the  joint 
criteria  of  minimizing  computation  time  and  restricting  the 
induced  error  in  the  solution  to  less  than  ten  percent. 

The  scheme  used  for  the  spatial  integration  was: 

1.  For  the  free-space  term,  and  therefore  the  perfectly 
conducting  term,  reduce  the  two-dimensional  integral 
to  a  one-dimensional  integral  of  a  finite  sum  of 
terms.  This  step  removes  the  difficulty  of  per¬ 
forming  numerical  integration  on  a  singular  inte¬ 
gral  without  introducing  an  error  which  is  large 
compared  to  machine  roundoff, 

2.  The  remaining  integration  was  performed  using  two- 
point  Gaussian  quadrature. 

3.  The  lossy  half-space  terms  were  integrated  by  a 
product  rule  of  Gaussian  quadrature,  i.e.,  a  three- 
point  rule  for  one  direction  of  integration  and  a 
two-point  rule  for  the  other  direction. 
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The  final  consideration  was  the  generation  of  the  lossy 
half-space  terms.  Their  generation . by  numerical  integration 
proved  to  be  prohibitively  expensive  thus  forcing  us  to  develop 
an  interpolation  scheme  for  these  terms.  We  found  that  even 
in  severe  tests  of  our  scheme,  a  30  by  30  interpolation  table 
proved  sufficient  to  obtain  the  same  final  answer,  to  within 
one  percent,  as  the  same  case  using  direct  calculation  of  these 
terms . 
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APPENDIX  A.  PROOF  THAT  F„  .  IS  INDEPENDENT 

N_ao 

s 

OF  THE  SCATTERING  OBJECT 


For  convenience  we  begin  our  proof  by  defining  the 
quantity 


AJl  1  s ,  2.-1  +  Is£,£-+1  1  2 


i  N-l 


(A-l) 


where  Ig£  k  is  defined  in  Equation  10.  The  definitions  for 
Z  =  1  and  Z  =  N  require  further  explanation.  For  Z  =  N, 

IgN  equals  the  limit  of  the  integral  over  a  hemispherical 
surface  having  a  radius  approaching  infinity.  Because  of  the 
radiation  condition  (or  exponential  damping  if  en  is  lossy) 
imposed  on  the  fields,  this  limit  is  zero.  The  case  Z  *  1  is 
treated  differently  for  the  horizontally  stratified  medium  and 
the  radially  stratified  medium.  For  the  horizontally  strati¬ 
fied  medium,  Igl  Q  equals  the  limit  of  the  integral  over  a 
hemispherical  surface  having  a  radius  approaching  infinity;  it 
is  zero  for  the  same  reason  as  IgN  N+1*  For  the  spherically 
stratified  case,  an  Ig^  ^  never  really  exists.  Using  this 
information,  we  complete  the  presentations  for  the  A^'s  for 
both  the  horizontally  stratified  and  the  radially  stratified 
medium  situations  as  follows; 


A 


1 


'si, 2 


(A- 2) 


=  IsN,N~l  (A_3) 

In  arriving  at  the  expressions  contained  in  Equations  A-2  and 
A-3,  it  was  necessary  to  use  the  boundary  conditions  at  infinity 
satisfied  by  the  fields.  In  order  to  obtain  the  following 
relationship  which  plays  an  important  role  in  this  proof 
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IVIUWIP 


k=  j+1 


(A- 4) 


Isj,k  “  _Isk, j 


it  will  be  necessary  to  employ  the  boundary  conditions  satis¬ 
fied  by  the  fields  at  each  interface.  Combining  Equations  8 
and  10,  we  have 


I  .  .  =  /  m . ,  dS 

sj,k  -  J  Tk 


where 


(A- 5) 


"j*  '  V|Sj^3jb-3jb«j.| 

Using  the  fact  that  for  general  vectors  A  and  B, 


ft 


•[AXB|  •  ft  |AtxBt  ] 


(A-6 ) 


(A- 7) 


where 


A t  =  A  •  ^-nnj 


(A- 8) 


we  can  rewrite  Equation  A-6  as 


*  fijk  ‘  [StjaxStjb-Et  jb^tja] 


(A- 9) 


Now  the  boundary  conditions  at  the  interface  are  such  that 


— tja  =  §tka  a=a'b 


(A-10) 


and 
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Using  the  fact  n .  ,  =  -n,  .  together  with  Equations  A-9,  A-10, 

j  r  &  K  /  j 

and  A-ll,  it  follows  that 


mjk  =  CA-12) 


This,  combined  with  the  fact  that  S.  .  and  S.  .  denote  the 

J  K  ,  j 

same  interfacing  proves  Equation  A-4.  Using  Equations  A-l 
through  A-4  and  forming  the  sum 


S  = 


(A-13) 


we  can  conclude  that  S=0,  or  equivalently  that 

N 


(A— 14) 


Since  J-  =  0  for  l  ?  N  ,  the  definition  of  A  in  conjunction 

X  cl  S  X» 

with  Equation  9  implies 


Ai  -  •Ivt*b  W, 


(A-15) 


Rewriting  Equation  15,  which  appeared  in  Section  2,  as 


and  substituting  Equations  A-14  and  A-15  into  this  equation, 
we  obtain 


N  ab 
s 


N 

-  £ 

£=1 


Vilab 


(A-17) 


In  order  to  draw  our  desired  conclusion,  it  is  only  necessary 
to  show  that  1^^^  is  independent  of  the  scattering  object. 

To  do  this,  we  write  the  definition  of  this  quantity  given  by 
Equation  12  in  expanded  form  as  follows: 


VZab 


*  J  5s<E-£a)-ilb<£ 


)  dV 


(A-18) 


Zb 


The  quantity  E^(r,ra)  is  the  electric  field,  observed  at  a 

point  reV^ ,  produced  by  the  current  distribution; 

a  =  a5(r-ra),  with  the  scattering  object  considered  to  be 
s 

absent;  and  izb  is  a  rigid  source  unaffected  by  the  scattering 
object;  Thus,  it  follows  that  *VJlab  is  independent  of  the 
scatterer,  and  from  Equation  A-17  it  follows  that  FN  ^  is 
independent  of  the  scatterer.  s 
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APPENDIX  B.  PROOF  THAT  Sj(r,r')  *  £k(r  Jr) 

The  proof  is  started  by  taking  the  curl  of  Equations  39 
and  47  and  then  substituting  Equations  40  and  41  into  the 
results  to  obtain 

V  x  (£»£.' )  =  iwu^g^rrr' )  (B-l) 

3 

V  x  (r,r* )  =  -iioe^  (r,r' )  +  S(r-r’)  ^2  ai^i 

=  -iuje^  (r  ,r' )  +  6  ( r- r ' )  J  (B-2) 

We  now  consider  the  special  cases  in  which  the  source  point 

r'  is  restricted  to  lie  in  the  volumes  V.  or  V.  ,  and  we  intro- 
„  3  k 

duce  the  notation 


r ’  =  r  eV. 

-a  ] 

(B-3) 

r '  =  r.ev, 

—  -b  k 

(B-4) 

Next  we  take  the  dot  product  from  the 

right  with  the  vector 

a  *  a  £=j 

(B-5) 

a  =  b  £=k 

(B-6) 

with  both  sides  of  Equations  B-l  and  B-2  to  obtain  Equations 
1  and  2  after  making  the  identifications 
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Combining  Equations  B-9  through  B-12,  it  follows 

XVjba  =  XVkab  (B-15) 

Using  Equations  5  and  12  as  well  as  Equations  B-5  through  B-9 
to  represent  Equation  B-15  in  more  detail,  and  evaluating  the 
volume  integrals  using  the  properties  of  the  5  function,  we 
obtain 


a*gj (£a,rb) • b  =  b*gk (rfa,ra) *a  (B-16) 

Using  general  properties  of  the  vector  dyadic  product,  we  can 
write 


a;Qj  (£a'Eb>  =  a-<*k(rb,ra) -b  (B-17) 

Using  the  fact  that  a  and  b  were  arbitrarily  chosen  vectors 
and  changing  notation  so  that  r&  ■  r  and  £b  =  r ' ,  we  obtain 
the  desired  result 

Gj (£,£')  -  2k(r',r)  reVj  r'eVk  (B-18) 
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APPENDIX  C.  OUTLINE  OF  ESSENTIAL 
DETAILS  TO  DERIVE  EQUATION  38 


An  essential  difference  in  obtaining  Equation  38  as  opposed 
to  Equation  37,  which  was  derived  in  detail  in  Section  2,  is  the 
choice  of  the  source  J^.  In  order  to  derive  Equation  38,  we 
choose 


1=1, 2, 3  (C-l) 


with  -aeVk,MNs;  as  opposed  to  the  choice  expressed  by  Equa¬ 
tion  16.  The  rigid  sources,  are  the  same  physical  sources 

for  the  excitation  as  before.  We  now  write  Equation  9  as 


Vkba  “  IVkab 

(C-2) 

-IVN  ab  ~  Is 

(C-3) 

hSL  Iv£ab 


W,2,-,N  (C-4) 


with  the  meaning  of  the  Aj's  defined  in  Appendix  A. 
in  that  appendix 

N 

S-S-o 

i-i 


As  seen 


(C-5) 


independent  of  our  choice  of  sources  J^a  and  J^*  Combining 
Equations  C-2,  C-3,  C-4,  and  C-5  we  obtain 
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N 

XVkba  Xs  "  2  IVjab 

j-1 


To  further  interpret  this  equation,  we  note  that 


XVkba 


-kb (^a} 


(C-6) 


(C-7) 


and  as  in  Section  2,  a  can  be  sequentially  chosen  as 
and  a3?  We  can  now  rewrite  Equation  C-6  as 


F.  w  +  I  =  a-E..  (r  ) 
kab  s  —kb  —a 


(C-8) 


kab 


N 

-  E 

j=l 


(C-9) 


Vjab 


and  Equation  C-8  can  be  related  to  Equation  14.  In  fact,  using 

the  evaluation  of  FN  ab  found  in  Appendix  A  and  comparing  that 

s 

expression  to  the  expression  given  for  Fkab  given  by  Equation 
C-9,  one  might  conclude  that  the  left-hand  side  of  Equation 
C-8  is  identical  to  the  left-hand  side  of  Equation  14.  This 
is  not  quite  the  case  because  the  implicit  dependence  on  r 

a 

causes  the  essential  difference.  In  Equation  14,  r  eV.T  , 

-a  Ns 

while  in  Equation  C-8,  r  eV.  .  Despite  this  difference,  the 

—a  ic 

same  arguments  that  were  used  in  Section  2  to  arrive  at  Equa¬ 
tion  31  from  Equation  14  can  be  used  to  arrive  at 


Ekb(ra>  "  ITkb(^a 


r  )  +  f  J(r)  • 


SN  (r^JdS 


raeVk  (C-10) 
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when  they  are  applied  to  Equation  C-8.  Using  Equation  B-15, 
we  can  write  Equation  C-10  as 


-kb(-a}  =  ^Tkb( 


r)  +  J  G,  (r  ,r)  *J(r)dS 


(C-ll) 


Recognizing  that  k  is  a  dummy  index,  which  for  completeness 
we  now  call  l,  and  making  the  same  notation  change  that  was 
used  in  Section  2  to  arrive  at  Equation  37  from  Equation  36, 
we  can  now  write  Equation  C-ll  as 

E£(r)  =  +  y*§£  <£'£')  *£(£'>  dS'  (C-12) 

*  s 

which  is  Equation  38 ,  and  thus  completes  this  appendix. 


